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Time domain step response expansions with coefficient algorithms are developed from the system 
function of a doubly terminated skin effect lossy coaxial transmission line. Three models of series 
impedance valid for (a) high frequency, (b) high and intermediate frequency, (c) high and low frequency, 
are incorporated into the system function. The system function is expanded via the method of polyno- 
mial expansions of analytic functions through generating relations and inverted term by term into the 
time domain through the inverse Laplace transform. Step responses for time domain reflectometry and 
transmission are developed, computed, and compared with experimental results. 
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1. Introduction 

In the analysis of skin effect lossy transmission lines one encounters propagation, characteristic 
impedance, and reflection functions which depend on Vs. For example, exp {—xy(s)} 9 Z (s), and 
p{s)= (Z — Zq(s))I(Z + Zo(s)) depend on terms like KV~s or R + KV^. Such functions appear 
in the study of pulse transmission in coaxial cables and integrated circuit transmission line struc- 
tures, etc., see references [1, 2, 3, 4, 5, 6]. 1 

A method of solution in the class of problems for R 9^0 is to affect rearrangements of the power 
series representing 7(5), Z (s), and p(s). A series representation for exp (-xy(s)) is realized 
through the technique of polynomial expansions via generating relations. Term by term Laplace 
transform inversion is then possible and the resulting functional series with certain expansion 
coefficient algorithms is the time domain solution. 

The time domain expansions which were obtained in this paper were specifically developed 
for the study of the pulse response of coaxial cables [6]. The conclusions of that study demonstrated 
that the planar skin effect transmission line series impedance models, the computational results, 
and the experimental data correlated quite well. 

We are now ready to pose the principal problem: determine the pulse distortion of a doubly 
terminated coaxial cable with skin effect losses. A circuit model of the system appears in figure 1 
for a transmission line of length / with 

(a) Zso(s), the source impedance; 

(b) Z r (s), the receiving end impedance; 

(c) Ei(s), the voltage step generator; 

(d) £"2(5, 0) and E 2 (s, /), step responses at input and output transmission line terminals 
respectively. 



""Electromagnetic Division, Pulse and Time Domain, NBS Boulder Laboratories, Boulder, Colorado 80302. 
1 Figures in brackets indicate the literature references at the end of this paper. 
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FIGURE 1. Circuit model of system function. 
Proceeding to the system function [6] 

E 2 (s,x) = Zq(s) 6 -^U) + Pr ( 5 ) e -(2/-^)^) 

Ei(s) Zo(s)+Z s0 (s) l-p s0 (s)pr(s)e- 21 ^ 

observe the voltage reflection function definitions 



(1.1) 



Pso(s) = 



Zso(s) —Zp{s) 

Z s0 (s) +Z Q (s) 



and 



Pr(s) 



Z r (s) —Zq(s) 

Z Z r (s)+Z (s) 



(1.2) 



Now y(s) and Zq(s), propagation and characteristic impedance functions respectively, define as 
follows [6]: 



y(s) = \/Z(s)Y(s), 



Z (s)- 



' Z(s) 
Y(s) 



(1.3) 
(1.4) 



The functions Z(s) and Y(s) which are unit length transmission line parameters are known as 
series impedance and shunt admittance respectively. They define via the relations [6] 



and 



Y(s)=Cs 
Z(s)=Ls+Ks ll2 +R 



(1.5) 
(1.6) 



where K denotes the planar skin effect parameter. 

The purpose of this paper is to obtain time domain step response expansions and expansion 
coefficient alogorithms from the system function (1.1) for three pairs of terminating impedances: 
(a) Z 80 (s)=0 and Z r (s)=Z (s), (b) Z 80 (s) = VLfC and Z r (s)=+oo, ( c ) Z s0 (s) = Z r (s) = VTJC. 
We are interested in step responses at x= and/ or x = l for small and large time arguments. 

In conjunction with the last remark observe that form (1.6) contains three models for series 
impedance [6] on selecting the proper value of/?: (a) R = 0, high frequency; (b) R = /?'/4, high and 
intermediate frequency; (c) R = /?', high and low frequency. R' is the direct current (dc) resistance 
of the inner and outer conductors per unit length. 

2. Properties and Expansions of y(s) and Zo(s) 

The following properties are important for extracting time domain expansions in later sections. 

(a) Singularity property. Let L x be the defined s domain of the inverse Laplace transform i.e., 
— 7r < arg 5 ^ 77. In view of (1.3), (1.4), (1.5) select the positive branch of s 1/2 . Now keeping in mind 
that L>0, X^O, R^O observe Ls + Ks 1/2 + R = possesses roots s* 12 and s] /2 in the left hand 
plane. Therefore, arg SkiLi(k = or 1) and y(s) (or Z (s) ) has one branch point 5 = 0. 
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(b) Bound on |e _/y(s) |. We want to show that | e /yU) | < 1 for all interior points in the region 
comprised of the right-half plane (rhp). According to Titchmarsh [7] the maximum- modulus theorem 
is valid for either branch of exp {— ly(s)} provided |exp {— ly(s)} | is single valued at the origin. 
To see this is true let 

y(8 + ./a>)= \y(8 \-JQ>)\€-M**°\ 

Then clearly exp (— |y(54 joj) | cos <£>(o\ co)) — > exp (0) = 1 independently of the path. Now on 
the rhp boundary \e ly{s) \ has a maximum equal to 1 and a mimimum equal to 0. Therefore, 
|e-W«)|<l for Re s>0. 

(c) Lower bound on ReZo(s). Let Z {) (s) = \Z () (s) \ e j % ,{s) and observe from (1.3), (1.4) and (b) 
above that — tt/2 ^ <p(8, co) ^ tt/2 implies — tt/2 ^ ^o(8, w) ^ ir/2. Hence, Re Z (§ + 7'co) = |Z (8 
-h 7"o>) I cos <po(o\ co) ^ 0. 

(d) Expansions. It is convenient to set slVLC=ar in obtaining a binomial expansion of 7(5). 
Then 

ffo + foo- 1 ^ *^ 
cr 



*y(*)-y(*)=<r|i+2 (- 1 )" " 



1 



(n+i)!' 



po=j-lVLC, p 2 = jl>iHLCyi\ 



Implementing the binomial theorem in (2.1) produces 

y(<r) = cr f 1 4- c o- 1 {^! ) + dftcr 1 ' 2 } - CiO-^dRJ + fl^cr 1 ' 
I 



-f 4;.]o- 2/2 } + 



+ (-i) n ™ 

. n + m + 1 



2 2, d n+m,k cr 



fc/2 



A- = 



(2.1) 



(2.2) 



(2.3) 



/here 



d ( n= — - 



«"*«- 



(/i~/c)!k! 

f K I K 2 R 1 2 
This series converges absolutely when » < |s| < S2 =M 2£ + y ^ + J J i.e., a proof of the bound 

appears in appendix A. Hence, rearrangement of the series (2.3) is realized by recognizing that 

n + 1 terms contribute cr ^ ( o" 11 ^ j. We require that integral value of A in (2.3) which cancels 

g-m (cr-m-i/2) and obtain the rearranged series 



with coefficient algorithms 



y(<r)=o-fl + |)a^ +1 o-¥}, 

m ■■ 
a[% +1 = (-D" J) (-l)» i c„ + ,4' + ) » 1 + 1 , 2 » +1 - 



(2.4) 

(2.5a) 
(2.5b) 
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(ID "I 1?2\ 
-y' TjFi ) 

i.e., valid for small |s|, yields the series 

y(s) = (RCyW*\l+ J (-l)«c„{|5^ + |s]" +1 }- (2.6) 

Proceeding as in (2.3) we find the rearrangement 

I , n+l 

y(s) = JRCM 1 + J <x<% + /*~), (2.7) 

with coefficients algorithms 2 

< ) 2n=(-l)"- 1 i (-l)*"c»+m-i4 2 .U«-m, (2.8a) 

m=0 

4 2 , ) 2 » + i= (-1)" 2 (-l) m c n+m ^ m+1;n _ m , (2.8b) 

m=0 

3. General Method for Time Domain Solutions of the System Function 

The following two inverse Laplace transforms are well known [8]: 

J?- 1 {-T^-} = ^) n,2 I n erfc [~~r-^J (zi = 0, 1, . . .), (3.1) 



Kty 






(3.2) 



Note k\ and k 2 are positive; I n erfc u(a — kj2t x ' 2 ) is the nth repeated integral of erfc # and possesses 
the recursion relation [9] 



.-/■ 

J u 



I n erfc a— I /*' ! erfc x dx 

J u 

where 1° erfc jt^erfc x. The notation H n {u) denotes Hermite polynomials of degree n. 

With (1.1), (3.1), (3.2) in mind we require a method of expanding A (w) exp (ag(w)) such that 
a z =a(x), and g(w) yields a Taylor expansion in cr~ 1/2 (or s l/2 ). From Boas and Buck [10] the 
technique of polynomial expansions of analytic functions through generating relations suffices in 
view of the following definition: 

A set of polynomials has a Sheffer representation if it is generated by the formal relation 

A(w) exp (ag(w)) = J Pn(a)w* (3.3) 



2 Superscript values of 1 and k s* 2 correspond to high and low frequency expansions respectively on all coefficients. 
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A(w) = J AnW«, /*„*(), 



g(w)=^ g„W", gl+0 . 



(3.4) 



If g\=0 and g„ =^ 0(n = 2, 3, . . .) we define P h {ol) to be a modified Sheffer polynomial. This 
definition clearly indicates Sheffer polynomials are rather general since they are functions of the 
expansion coefficients in A(w), g(w), exp (u). The mechanics of obtaining P n (ot) are taken up 
in section 4. 

The application of (3.3) to (1.1) requires a hierarchy of coefficient algorithms such that recursive 
properties hold within and between algorithms. To clarify this point note the d coefficients of (2.5) 
require recursivity within (2.5) to determine the "a" coefficients while application of (3.3) to exp 
(— y(cr)) yields at least recursivity between a new algorithm and (2.5). 

4. Ideal Step Response 

A very simple pair of impedance values for (1.1), Z } (s) =Zo(s) and Z s o(s) =0, leads to the 
o~ domain ideal step response 

E 2 (aJ)=E^^~, E = ElV7I;. (4.1) 

a 

By virtue of (3.3) this expression becomes 

&(<r, /) -&-% CXP (-<°- +a ^°- l/8 > f ,(<r, I), (4.2) 



f,(<r, J) = expf+ J rf*r~W|, 



(4.3) 



where g\}l^ = — a{)l (n ^ 3) and (4.3) is a member of (3.3). According to Knopp|ll|the main 
rearrangement theorem applies to (4.3) since 

i.e., (2.4) is absolutely convergent. Hence, proceeding to a tractable expansion of (4.3) obtains the 
rearrangement 






V-(x)]' = 



n = ' m = 



}*-<r-W J G<i> B>1 o-»/' (4.4) 



where all G coefficients 3 define [12] through the recursion 



G V», i = ZZT)2{P( k + 1 )- »>«W.-* ■ (4.5a) 



3 The third subscript of G means its integer value is a multiplier of each g coefficient on the right-hand side of the definition equation (4.5a). It corresponds to the 
number of reflections (or transmissions) which the input step voltage experiences. 
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with initial values 4 

G&i-l. Cg>..,- (*?>)*. Ci%., = «i 1} . (4.5b) 

From (4.3), (4.4), and (4.5) we obtain a double power series which convolutes with respect to 
subscripts and yields 

f,( ff ,!)=fif>"*- (4.6) 

H = 

n (7(1) 
1,(1) = V ^fc,n-fc,l L(i) _. i 5 

such that all coefficients b\}\ are Sheffer polynomials with argument x = I. 

Now following Doetsch [13] term by term inversion of (4.2) is permissible. Hence, in view of 
(3.1), (4.2), and (4.6) the time domain expression becomes 

e 2 ( T , /) = Ee~*% J 6« { 4 ( T - l)}'" 2 /» erfc {^ (t- 1)-*]. (48) 

n = ^2 J 

X 

To determine the t interval of uniform convergence observe that V |6^\ |S^" n ^ 2 (S4 > S 2 / VLC) 

w=0 

is a convergent series. If we require term by term dominance of the expansion in (4.8) the inequality 



{ 4(r _ l)}n/2/„ erfc feii ( T -_ !)-l/«| ^ S n/2 



(4.9) 



determines the interval 1 ^ r ^ To where To is that value of r which makes (4.9) an equality. 

To compute the step response for large time arguments revert to (2.7) and consider the s 
domain expansion with coefficients g$ = — lVRCa[ 2 ^ m+1 (m = 0, 1 . . .). Then 



IVRCs" 



E t (s, l)=E - = E exp 2 ggs 2 . 

5 S ^m=0 > 



(4.10) 



For < \s\ ^ So we use the argument supporting (4.6) to obtain 

-iVRCs 1 '* x 

E t (s,l)=E- JftS?!*"*. (4.11) 

[n/2] C(2) 

where [n/2] denotes the integer function. The coefficients fc^i are modified Sheffer polynomials 
with argument x = L 

Now inversion of (4.11) through use of (4.3) and (4.4) carries \s\ beyond S . Hence, according 
to Doetsch [14], we can expect the time domain expansion to be asymptotic. Indeed, this turns 
out to be true under rearrangement of (4.11) into odd and even numbered terms from a proof in 
appendix B. Hence, the conclusion contains the form 



4 Initial values of G coefficients assume the same form as (4.5b) when superscripts are changed. 

5 We define b[ k) n ■ 1, independent of n and k. 
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e 2 ( t ,/)~£erfc(-^ r )+^ i ex„ -— X(2*g, 2 L, n/2 (4.13) 



M = l 



where the operator ( ) denotes rearrangement with respect to n. 

5. Open Circuit Step Response 

To handle the open circuit case revert to (1.1), set Z s q(s) = y/LjC and Z r (s) =-f oo, and observe 
the input terminal step response 

At this point it is worthwhile to investigate the possibility of expressing the r time domain solution 
of (5.1) as a sequence of reflections. This corresponds to expressing (1— p s{ )(s)e~ 2lyis) )~ l as a 
geometric series. The condition \p s o(s)e~ 2ly{s) \< 1 is satisfied in view of 3b and 3c. We remark 
that it is also quite easy to calculate an upper bound for |p. s -o(s) | by Stannard's method [IS]. 
Applying the geometric expansion to (5.1) yields 

Zo(or) + VL/C ^ (O.-i) 

Following the expansion technique for y(cr) obtains the "unreflected term" representation 

*».•<"• °> B > t \°+lmr = l + i ^^» g " ( ^ 1) » (5 - 3a) 

Zo(o-) + VL/C 2 ~ 
<#> 2B =(-l)»-« 2 (-l)*c»+*<ftV*.«, (5.3b) 

A = 

ai\\n +i = (" 1) " 2 (- 1 ) A '<w + . <*&*+ . , 2A-+ ■ • (5.3c) 

/f = 

Integral powers of voltage reflection functions expand into the form 

» / m+n \ 

p« (o-) = (-2)« £ Ailli,-i<r~{~' (n=h2, . . .) (5.4) 

m=0 

which defines itself according to the alternate recursion relations 6 

m 

A(\) =a {l) and 4 (l) = V n {1) A {1) fo <%\ 

A- = 

We expand the exponent functions of (5.2) into the form 

e -**<*> =€" 2 KV a exp |-2p{o-f ajftcr 1 '*} ] j? b£] tp o- mf2 (5.6) 

1 Of course relations (4.5) and (5.5) give the same answer. However, their implementation points out one difference which we shall mention in section 7. 
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such that p = n, n+1. Hence, from (5.3), (5.4), and (5.6) the nth term becomes 

U 2 \l l bHl k , 2p Ai% n _ 1 }<r-{ !! T L ) 

m = A fc=0 J 



■f 2 ^.^Km)^^)) (»=o,i, 



•)• (5.7) 



The prime notation on sigma indicates the series (5.3a) replaces (5.7) at p= and, at p= 1, 6£|> 2 re- 
places the sum ^ &m-*,2^t+i, -i- 

/c = 

Hence, the time domain response yields 
e 2 (r,0)=£ | (-2)4 2'«- ?PB( 4i 2 { 2 ^Ut^Wi..-!} 

ffl + n K ^ m >, w + n + 1 

x{4(t-2/>)} 2 /«+» + £ 2 *&*.*4&.« {4(r-2p} • 
o o m = Jjfe=o |o 

X /m + n+l| er f c { pa (i) 1 ( T -2p)- 1 / 2 } | (5.8) 

1 2 

and the term for p = inverts to the series 

St..(r, 0) =-+ J aaU, /m ,o N - (5.9) 



A '• m+l r (" 



This series converges uniformly when 



(T*+>/*)/r('£fi'us« ^ (510) 



Clearly inequality (4.9) (n replaced by m + n) extends to cover (5.8). 
To determine the time domain asymptotic expansion let 

Zo(s) 



■y { l + wM+w-nlW { 1+ ± a[% +1 s<» - ')/2 j 1 1 + ^ 1/2 • (5.11) 

Now define a set of coefficients such that 

4 2 >=1, 4 2 >=-n --n 1 ' 2 , 4 2 »=T, 2 -r,i'%< 2 ; ) 1 , 4 2) =-t?| /2 o<,Vi (« = 3,4, ...). 
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Then (5.11) assumes the form 



Zo(s) 



Zo(s) + 



L I f Vos 112 £J 






c 



In a very similar manner we obtain the expansion 

1 



P80\S) 



- V A i3) 5 n/2 



fl/2 



with coefficient definitions 

4%=1- ^ n = r, n -2r,}J\ Ai% = 2{r, 2 -r,^af\}, 4 3 , ) o=-2^M%_ 1 
For convenience redefining the coefficients of (2.7) yields 



The exponential function of (5.1) has the expansion 



c A-/2 



e _ 2 , 7 ,, )= ^ 2 , x ^ 6Y,)„,,,»^= 2 H 3 > B/ 



with the coefficient algorithm 



ftff.- 2 



« /-;<3) 

r> °*,n-fc,2 



4! 



(i>>2). 



(5.12) 



(5.13) 



(n=3,4 . . .). 



af\ = -lVRC, «l"„ = -/VfiC<» )M (n=2,3, . . .). a<«, +1 =g<?> (n = 0, 1, . . .). 



(5.14) 



Combining expansions (5.13) and (5.14) obtains 

l-p„( S )e-"*'>= * 2 £ «?),-/«, (5.15) 

subject to coefficient algorithms 

M<, 2) = 2 , M< 2 > = tjo +i<i»> a + b\% , 

M ( » 2) = £ ^*,oH 3) 2 ("=2,3, . . .). 

A- = 

Expansions (5.12) through (5.15) constitute a scheme for the series solution of (5.1) by means of the 
two quotients 



E 2 (s,0)=E l (s) 



2 4?>** /2 

£ M,?V" 2 



jc n 



-f 6 -2/VKCs 12 



2 1 4W& ^ 

M =0()^ k=0 



J M»s n ' ! 



(5.16) 
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388-513 O - 70 - 3 



Since the power series above are analytic for < \s\ ^ So we may perform division termwise 
to obtain 



E 2 (s,0)=E 1 (s) 

such that recursion formulas 



,,=0 n = J 



(5.17) 



n £j n + k k 



WW* 



and 



W\? 



n-l 

2 

k = () 



^^l k b$ 2 -^M<Zl k W^ 



(n=l,2, . . .), 



with initial values JV\{ )= W { ^ ] — - enable all coefficient determinations. For a time domain expansion 

of (5.17) invoke the gamma function equation V (z)T (1 — z) = ttcscwz (z complex and nonintegral 
real) and observe 



2 



IVRC 



l-RC 



+ o erfc l^r +£e 4 ' 2 r » 



H n -\ 



(3). 



IVRC 

21 1 ' 2 



2"-Wnt" 



(5.18) 



as t — > oc is a complete asymptotic expansion, 
The short circuit step response with imp 
from the basic a and s domain expansions of this section 



The short circuit step response with impedance values Z r {s) = and Z s o(s) = \jy; is obtained 



6. Time Domain Reflectometry and Transmission 



We select the impedance values Z s o(s)= z Z r (s) = y-7; and require a solution of (1.1) for 
^ x ^ / and ^ t ^ T < oo. Let p r (cr) = p s o(<r) = p(cr) and define three relations 



a 2n (x) =2ln + x, 

ao n + x (x)=2l(n-\-l) —x, 

& n (x) =~j OL n (x), (/l = 0, 1, 2 . . .) 



which include two integer functions 



&n(0)=2 



n-f 1 

2 



and a n (l) =2 



+ 1 



at input and output terminals respectively. 

Invoking the geometric expansion in (4.2) via (6.1) yields 



(6.1) 
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7 I \ 

E 2 (o-..v) = E,(o-) - % p»(o-)€-«"W*W. (6.2) 

Zo(a) + > /§" =0 

To implement the rearranged series set 

T n (a, x) = exp I & n (x) j? gti^ " ("V^l (6.3) 

Let the definition equation 

w=o J m=o (6.4) 

analogous to (4.4) apply such that all 7 coefficients obey the relations (4.5). Then as in (5.5) we have 
the expansion 

r„(cr,x)=£ B\» n {x)<T"" 12 (6.5a) 

iw=0 

contingent on the definition 

B^, n (x) = 2 ak n(x)A%- k . (6.5!,) 

The convergence of (6.5a) is uniform over ^ x ^ / and °° ^ |cr| ^ S4. Since the inverse Laplace 
transform of (6.2) is zero when r= [7 T // VLC] + 1 (integer function) we do not inquire about the 1 
convergence of (6.2) with respect to //. 

For ^ x ^ / and ^ n ^ °°, the or domain step response of (6.2) with support from (5.3), (5.4) 
and (6.5) becomes 

E i (*,x)=E l (<r) %'I(-2)H "•<""'" ! { £ f| 2 4 1 ].- -,*#.(*)] 

w = , A=<) Jj m=0 \ } 

-(!LL!L\ x f A 1 _A + /, + i, 

^ - '+2 S4i ffl A(4 2 

f- — ft *"■**» — ft ■* ■* 



fc=0jm=0 Q 

Xexp {-&„(*) ((T+ after 1 ' 2 )} (6.6) 



where the polynomials of degree /: are Sheffer and the prime on sigma indicates two exceptions: 

(a) Replace £ ^.-...-.^.W with B#,(*) at «=0, 

m = 

(b) Replace (6.6) at n = and x = with (5.3a). 
Hence, the time domain response yields 

e 2 (r, x)=E §'/(-2)»« "***» f £ f| E ^i-m..-i«a. ) .(*) 
x{4(t-«„(*))}~/* + »+ J I V 4&_, B5i>.(*) 

X{4(T-a„U))} _ ^/ A ' l " ,, }ericj^4^ a V , 1 {T-a„U)}- |/ 4 (6.7) 

J I ^ J 



1 1 
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where (5.10) replaces (6.7) at n = and x = Q. 

To determine the T interval of uniform convergence invoke the inequality 



k + p 



{4(r -«»(*))} 2 / A+ "erfc 



^aft (t -&(*))-»} <S 



A- + j> 
4 2 



(6.8) 



for p = n, n+ 1. 

For large time intervals we first note that (5.13) and (5.14) furnish a starting point for the 
expansion 

-1/2 x 

1 - p»(*)e-«W> = , „ £ M<f>5»' 2 (6.9) 

U -r no5 ; w=0 

such that Mf ] coefficients define according to the rules: 



fc-0 

From (5-12) and (6.9) we have the product 



Z«(s) 



I 



(1 + T>0S 1/2 ) J ^V 1 ' 2 



/ 7 l _ Q 2(e) € -2ly(s) 



(6.10) 



Following the methods of (6.3) and (6.4) produces 

fc=0 <-f»=0 



6 -ap(x)Vi?CsV 2 j? B( n 2> p (^)5»/ 2 , 



(6.11a) 



contingent on the following definitions and properties: 

[n/2] /v/c(y\ 

B&<*) - 1, SS 2 >„« - 0, £<„%(/) - «&(/). 
We define the equations 



CJJW*) = 2 ^-.Wgoto (* > 0), 



(4« 



C(f,) (0)=^ 
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(6.11b) 
(6.11c) 



(6.12a) 



(6.12b) 



c< 2 >(*) = 2 2 ^8MfcUt*IBi(*). 



(6.12c) 



Q!i(*)-2 W„,,,(*)M</' 



Af„,A-(x)=- 



Ml« 



n> 1 
x^O 



(6.12d) 



From (6.11) and (6.12) three special cases condense into the forms 

/Vo.n(%)=/V ,,(*) = l/M< 3 >, 



Nt,o(x)=- 



M 3 > 



W».,(*) — 



^«+4fi-Af\ 3 VM , 3) 



«i» 



(6.13) 



%&0. 



In the s domain we have the representation 

£ 2 (5, x) = 2 (- 1 )*«~°*< ■ r)V ^ s,/2 2 ^».fc(*)* (B " 3)/2 



(6.14) 



which converges uniformly in x ands overO =£ a: =s / andO < 60 * \s | * So. It is necessary to consider 
three positions on the line: (a) x = 0, (b) * = /, (c) < x < I. At x = the terms corresponding to 
n = on Laplace transform inversion contribute t" 2 , a term which cancels on using the series 
expansion of/ 1 erfc ( ). Hence, observe the complete asymptotic expansion for time domain 
reflectometry (TDR). 



e 2 (t, 0) ~ E 



2NolVRC £ 



(-)' 



21VrC }'"-' 
t^ J 



'»=' m\T 1 + 



1 — m 



+ /v I)0 (o)f/(o 



/Vkc 



+ X (-l)"-W 2n ,o(0)r(^)/: v - '-JV^^erfc -„ 

n=l I * 

_mc x tl n -l\ I/2 i 

-e ' \2 /V»,,(0) 



(6.15) 



2"-2f(n-D/2 

as ( ^ °°. In view of (6.12) and the coefficient definitions of section 5 (6.15) at t = °° reduces to 



e 2 K0) 



» 



(6.16) 



\/^ + IR 



the d-c voltage divider equation. 

Over < x < I we use the expansion of/ 1 erfc ( ) again to cancel the term associated with 
t 1 ' 2 and extract the complete asymptotic expansion 
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e->(t, x)~E\ N VRC 2 (~ 1) A «*W 2 7 pr^T ^7T72 — 

*=o m=i mir 1+- 



2 



+i(-i)w».*«fe( fl * ( ^ c ] 



/7 _ | «A-U) VKC ' 



as f — » oo. In the case of time domain transmission (TDT) note %=/ means ao(l) = <X\ (/) and, 
consequently the coefficient of iVo vanishes. Also at t — °° and i = /we obtain the d-c voltage divider 
expression 

L 



r 
e 2 K/) = — = (6.18) 



7. Realization of TDR and TDT Data 



7.1. Theoretical and Experimental Results 

A transmission line with the following parameter values is illustrated: L == 1.184 X 10~ n H/m, 
C-2.46X10- 7 F/m, /= 137.2 m (450 ft), R' = 1.466 X 10" 2 fi/m, KVs~= 3.35 X 10" 5 V$~ H/m. 
Three curves of each figure correspond to three models of Z (s) (sec. 1). Each theoretical curve 
approximates experimental data best over a certain time interval and the relative error between 
the two sets of data (taking either set as correct) over that time interval is realized through the 
following example: 

TDT Data 

t (microseconds) 1 2 3 41 _, . R „„ /— 

% error 2 1 1 2J 4 

In view of a 3 percent uncertainty of the experimental system [6], a maximal 2 percent error for 
the appropriate Z(s) in the data of the figures below is optimal. Experimental data points for t less 
than 10 ns are not reliable [6]. 



7.2. Computational Observations 

Computer programs were written on a CDC 3800 digital computer in Fortran language to 
realize theoretical TDR and TDT data. We mention certain considerations of the computational 
procedure for (6.7) by starting with c and d il) coefficients of (2.5a). Their product possesses 
recursions 

»~5 (" + £--] (n-k) 

ITU) = -flj(l) /TO) =Ji AZ £2*T(l) 

*»+i,o n+l P»Tn,o, r n>k+1 (2i+i)(2Jfc+2) /So n ' k (7.2.1) 



vhere 



F^=c„ +k -,4^ fe)2fc (0 *£ k *£ n and n= 1, 2, . . .). 
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Similarly from (2.5b) 



1 



n(D =q -7)(i) 



7)(1) 



n + k + l) (n-k) 

zl & n(i ) 

(2A + 2)(2* + 3) /So n ' A ' 



(7.2.2) 



where 



D{ n\ = c n +kd { 1 ll k + 1 , 2 /c+ 1 (° ^ A: ^ n and ti = 0, 1, 



.). 



In computing B$ n {x) (at jc = 0, /) it makes no difference whether (4.5) or (5.5) is implemented 
to realize ./[J^^ since indexing occurs over k. However, suppose inquiring only about the magnitude 
of reflection n occurs. Then coefficients A^ n compute via (4.5) by recurring over the mth now of 
the nth column 7 whereas (5.5) requires recursion between columns n and n — 1. 

For the time interval 8 ^ t ^ 10 ^ts it is necessary to set k = 1(1)25 and n= 1(1)10 to compute 

the Sheffer polynomials 

k 

V A{\) R(D (r) 

m = 

Their ranges of magnitude at x = and x~ I appear below: 

k=\ /c-25 



n = 
n=\0 



/oXlO" 2 


f x x 10- 14 


/2X10- 16 


h x iO- 23 



K|/,|<6(i=0, . . ., 3). 



Now determining numerical values of erfc u 



7 (D 



(u^~jr 1 a n (x)X(r-~a n (x))-^ 



through its power series and using the recursion formula [16]. 

I p erfc u — I p ~ l erfc u + — I p ~ 2 erfc u 

P 2p 

enables term by term computation of (6.7). In reference to (6.7) let 

e 2 (r, x) =2 ^2,n(r, x) and e 2 ,n(r, x) =2 ^,«,a-(t, x) . 

Then a pair of practical terminating conditions for computing (6.7) can be stipulated: 



(7.2.5) 



(a) I e 2 ,m,q(m)(T, x) | = 10" 

(b) \e 2tP (T,x) | = IO 4 , 



m + q{m)^25, m=l(l)p; 
10. 



To clarify these conditions suppose that e 2 (r, 0) is wanted for r= 14.1 (t— 10 juls) . Reference to 

(6.7) shows that n— 1(1)14 while numerical computation reveals e2,m,g(m)(T, 0) satisfies condition 

(a) for ra^ 1 (1)4. At p = 4 conditions (a) and (b) are met and all computations on (6.7) are terminated. 

A practical value for the interval of convergence may be determined by invoking the inequality 



7 Note the subscript position reversal of A\^ n in comparing (> [ J \ ti of (4.5). 

H Setting k— 1 (1)50 and /?= 1(1)10 jumps computing time on the CDC 3800 from 30 s to 3 min. This happens because so many convolutions occur. However, 
it is not necessary to go beyond k ■ = 25 for this particular transmission line. 
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2*(T-a„ (*))*" /'eric (^ra H (x)(T-a n (x))- 1 l i ) ^2><t- 



l2 + V (7.2.6) 

Substituting the asymptotic expression \(z) ~ 2tt€ ~ z t z ~ 1/2 , A:=25, and computing that value 
of t which satisfies the equation 

T={r(25/2-hl)} 2 / 25 S 4 - 1 (7.2.7) 

produces t==7AXlO~ H seconds. Now £2, i,«(i)(t, x) |x=o, 1 does not satisfy condition (a) at 
£ = 2.1 X 10~ 5 s. Hence, (7.2.7) gives a pessimistic (but useful) guideline for the interval of conver- 
gence. 

In asymptotic calculations terms corresponding to No,k(x) and Ni,k(x) are dominant and no 
more then three terms of each series (odd and even) are needed for this combination of parameter 
values. Asymptotic results are used for £ = 20, 30, 40 /xs, 

8. Appendix A. Bound for |(j| 12 

We have as a condition from the binomial series (2.1) 

liSo + jSzor 1 / 2 ! < |cr|- (Al) 

Let |o~| = v 2 and invoke the triangular inequality to obtain 9 

i; 2 -|/3o-f/32a- 1/2 | 2*v 2 ~l32V-po=f(v). (A2) 

Calculating the positive zero of f(v) yields 



j3 2 + Vj8i+4flo 

vo = (A3) 

It remains to show that f(v) ^0 for v^Vq. Noting f(v) has one extremum, a minimum, at 

v = ~r and f (v) ^ for v ^ /W2 completes the argument. 

This bound is sharp since the selection or ll2 = vo is an admissible value. Transforming from 
(j to s produces the result 



ii + v^ + !F- < a *> 



9. Appendix B. Asymptotic Expansion Discussion 

The problem is to show that 



e,(t, /) ~£erfc (8t~^) +-|" 2 exp I -8*/"' ] f J fc£> + , 

77 I J I A T =0 



v H 2n (8f-"») H iH +1 (8t-M) . 

22»tV z ' 



9 The author is indebted to O. N. Strand for this argument. 
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is completely asymptotic ( 8 = - vRC I. 

Let (p n (t) = H u -i(8t~ h2 )lt nl2 and define a sequence {i//»(0}/i such that 

i. i/f/i+i(0 n r 

lim ~ — -— =0 tor every rc. 



(B2) 



Then (B2) forms an asymptotic sequence [17] uniformly in n. Using the definition of Hermite 
polynomials [18] observe the limit quotient 



lim 







2»+l 


<P2n+l(t) 
(p2n(t) 


~lim 

t-*ao 


t 2 

2h + 1 

t 2 



2 = 2 

6~S 



as t—> °°. Hence, a rearrangement is necessary to satisfy (B2). 

It is quite easy to see that odd and even subsequences satisfy respectively 



lim 



(P2jl+s(t) 

<P2n+l(t) 



lim 



^^UlimO(-) 



A complete asymptotic expansion with respect to {ifjk(t)}k obeys the relation [19] 



t— *°°. 



(B3) 



(B4) 



(B5) 



By definition of {<p 2 »(*)}2n, |?2jv(0l ^ |<p2n(*)|(n = iV, TV + 1, . . .) as £->°c. Also b$ +lfl origi- 
originates from a convergent series, i.e.. b 2 2 n+i i ^ A?< °°. We have 



V />(2) ?*»(*) 

2^ "2H+1, l n2n 



^Af 



^2.v(*) 



1 



2^ = ° (^v(O). 



(B6) 



as t — > oo, a relation which satisfies (B5). Clearly (6) also holds for the odd numbered terms. This 
argument extends to (6.15) and (6.17) in view of the fact N n , k(x) is uniformly bounded. 



The author is sincerely grateful to Prof. N. S. Nahman for his guidance and inspiration through- 
out the course of this project. 
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